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Multi-scale  and  Multi-physics  Numerical  Methods  For  Modeling  Transport  in  Mesoscopic  Systems  (a  proposal 
submitted  to  Numerical  Analysis  Program,  Mathematical  Sciences) 

ABSTRACT 

In  this  project,  we  have  accomplished  in  development  of  algorithms  to  model  transport  and  eletromagnetic  processes  in  mesoscopic  systems 
such  as  nano-electronics  and  biological  membrane,  and  layered  inhomogeneous  media. 

Specifically,  the  following  results  have  been  obtained  resulting  in  the  publication  of  6  peer-referred  journal  papers  and  a  third  part  of  a 
Cambridge  University  Press  book. 

(1)  fast  integral  solver  for  quantum  dots  in  3-D  layered  media.  The  fast  solver  is  based  on  a  window  accelerated  method  for  computing  the 
layered  Green's  function  and  wide  band  Fast  multipole  methods  for  Hankel  waves. 

(2)  a  new  linear  scaling  discontinuous  Galerkin  density  functional  theory,  which  provide  a  brand  new  approach  in  combining  physics-based 
orbitals 

and  piece-wise  polynomial  finite  element  basis  in  finding  the  ground  state  energy  of  the  DFT  for  quantum  systems. 

(3)  numerical  methods  for  computation  of  electrostatics  in  ion-channel  transport, 

(4)  a  new  parallel  solver  for  elliptic  PDEs  by  combining  random  walk  Feynmann-Kac  formula  and  local  boundary  integral  equations  for 
extreme  computing, 

(5)  an  improved  device  adaptive  inflow  boundary  condition  for  Wigner  quantum  transport  equations. 

Also,  a  book  titled  "Computational  Methods  for  Electromagnetic  Phenomena:  electrostatics  in  solvation,  scattering  and  electron  transport" 
was  published  by  Cambridge  University  Press  on  Feb.  25,  2013.  The  work  on  electron  transport  (Part  III  of  the  book)  results  from  this 
project. 
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Scientific  Progress 

[1]  On  fast  integral  solver  for  quantum  dots,  we  have  developed  a  parallel  fast  algorithm  for  computing  the  product  of  the 
discretized  Helmholtz  integral  operator  in  layered  media  and  a  vector  in  0(NqNzA2  Nx  Ny  logNx  Ny)  operations.  Here  Nx  Ny  Nz 
is  the  number  of  sources  and  Nq  is  the  number  of  quadrature  points  used  in  the  evaluation  of  the  Sommerfeld  integral  in  the 
definition  of  layered  Green's  function.  The  fast  solver  is  based  on  two  important 

techniques  which  reduce  the  cost  of  quadrature  summations  in  the  Sommerfeld  contour  integral  for  Green's  functions  in  3-D 
layered  media. 

The  first  technique  is  the  removal  of  surface  pole  effects  along  the  real  axis  integration  contour  by  identifying  the  pole  locations 
with  a  Discrete  Wavelet  Transform;  In  the  second  technique,  we  apply  a  window-based  high  frequency  filter  to  shorten  the 
contour  length.  As  a  result,  the  integral  operator  for  the  3-D  layered  media  can  be  efficiently  written  as  a  sum  of  2-D  Hankel 
cylindrical  integral  operators,  and  the  latter  can  be  calculated  by  either  a  tree-code  or  a  2-D  wideband  FMM  in  a  fast  manner. 


[2]  On  the  linear  scaling  quantum  DFT  algorithm,  we  have  proposed  a  new  framework  using  discontinuous  Gakerlin  method  for 
linear  scaling  methods  for  density  functional  theory,  which  forms  the  fundamental  approach  for  studying  quantum  system 
ground  state  energy.  The  salient  feature  of  this  framework  is  the  flexibility  of  using  hybrid  physics-based  local  orbitals  and 
accuracy-guaranteed  piecewise  polynomial  basis  in  representing  the  Hamiltonian  of  the  many  body  system.  Such  a  flexibility  is 
made  possible  by  using  the  discontinuous  Galerkin  method  to  approximate  the  Hamiltonian  matrix  elements  with  proper 
constructions  of  numerical  DG  fluxes  at  the  finite  element  interfaces. 


[3]  Computation  of  electrostatics  in  ion-channel  transport. 

Here  we  investigate  the  numerical  accuracy  and  efficiency  in  computing  the 
electrostatic  potential  for  an  ion  channel  model  made  of  finite-height  cylinder 

embedded  in  a  layered  dielectric/electrolyte  medium  representing  a  biological  membrane  and  ionic  solvents.  Two  numerical 
techniques,  a  specially  designed  boundary  integral  equation  method  and  an  image  charge  method,  have  been  investigated  and 
compared  in  terms  of  accuracy  and  efficiency  for  computing  the  electrostatic  potential. 

[4]  A  new  parallel  solver  for  elliptic  PDEs  for  extreme  computing. 

In  this  result,  a  hybrid  approach  for  solving  the  Laplace  equation  in  general  3-D  domains  is  proposed.  The  approach  is  based  on 
a  local  method  for  the  Dirichletto-Neumann  (DtN)  mapping  of  a  Laplace  equation  by  combining  a  deterministic  (local)  boundary 
integral  equation  (BIE)  method  and  the  probabilistic  Feynman-Kac  formula  for  solutions  of  elliptic  partial  differential  equations. 
This  hybridization  produces  a  parallel  algorithm  where  the  bulk  of  the  computation  has  no  need  for  data  communication 
between  processors,  therefore  it  has  great  potential  for  highly  scalable  solver  for  elliptic  PDEs  for  extreme  scale  computing. 


[5]  An  improved  device  adaptive  inflow  boundary  condition  for  Wigner  quantum  transport  equations. 

An  improved  inflow  boundary  condition  is  proposed  for  Wigner  equations 
in  simulating  a  resonant  tunneling  diode  (RTD),  which  takes  into  consideration 

the  band  structure  of  the  device  and  the  effect  of  the  quantum  interaction  inside  the  quantum  device.  Numerical  results  on 
computing  the  electron  density  inside  the  RTD  under  various  incident  waves  and  non-zero  bias  conditions  show  much 
improvement  by  the  new  boundary  condition  over  the  traditional  Frensley  inflow  boundary  condition. 
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